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Abstract 

This  paper  describes  measurements  and  numerical  simulations  of  the  movement  of  longshore  sand  waves  observed 
at  Southampton  Beach,  Long  Island,  New  York,  These  large  morphologic  features  are  created  by  periodic  opening 
of  a  small  inlet  and  subsequent  welding  of  its  ebb  shoal  to  the  beach.  Longshore  sand  waves  are  wave-like  forms  that 
maintain  identity  while  moving  along  the  shore  and  represent  a  simple  type  of  collective  sediment  movement  or  large- 
scale  coastal  behavior.  Despite  being  described  as  early  as  1939,  few  measurements  of  longshore  sand  waves  have 
been  made,  and  their  dynamic  behavior  is  not  well  understood.  Eleven  longshore  sand  waves  were  identified  along 
Southampton  Beach  and  have  an  average  length  of  0.75  km  and  an  amplitude  of  40  m.  Migration  speed  was  greatest 
under  winter  wave  action  and  reached  an  average  1.09  km  year"*  and  yearly  average  speed  including  summer  calms 
was  a  much  lower  0.35  km  year”^.  During  a  16-month  observation  period,  the  longshore  sand  waves  did  not  disperse 
during  their  steady  migration  westward  in  the  direction  of  predominant  wave  incidence,  and  longer  waves  moved 
faster  than  those  with  smaller  wavelengths.  A  numerical  model  of  shoreline  change  including  both  standard  particulate 
longshore  sand  transport  and  a  component  describing  migration  of  a  longshore  sand  wave  is  developed  in  which  the 
velocity  of  the  wave  is  related  to  the  calculated  longshore  discharge  of  water.  The  model  shows  wave  refraction  to  be 
an  important  mechanism  contributing  to  the  longevity  of  longshore  sand  waves,  acting  to  retard  the  otherwise 
expected  diffusion.  Model  calculations  of  the  migration  of  a  single  longshore  sand  wave  agree  quantitatively  with  the 
observations. 


1.  Introduction 

This  paper  describes  measurements  and  numeri¬ 
cal  simulations  of  the  movement  of  longshore  sand 
waves  observed  at  Southampton  Beach,  on  the 
south  shore  of  Long  Island,  New  York.  Longshore 
sand  waves  are  wave-like  forms  that  maintain 
identity  while  moving  along  the  shore.  We  associ¬ 
ate  the  word  “longshore”  with  these  sand  bodies 
to  distinguish  them  from  other  morphologic  wave 
forms  that  appear  in  navigation  channels,  rivers. 


and  coastal  waters,  but  which  are  not  directly 
related  to  shore  processes.  Longshore  sand  waves 
are  relatively  large  features  and  have  been  docu¬ 
mented  to  persist  over  decades  (Bakker,  1968)  and 
even  over  a  century  (Verhagen,  1989),  as  opposed 
to  beach  cusps,  which  are  smaller  features  with 
much  shorter  lifetimes  on  the  scale  of  changes  in 
wave  conditions.  Other  terminology  has  been  used, 
such  as  “cusp-type  sand  waves”  (Sonu,  1968), 
“migrating  sand  humps”  (Bruun,  1954;  Grove 
et  al.,  1987),  “accretion  and  erosion  waves” 
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(Inman,  1987),  and  “sand  waves”  (Verhagen, 
1989).  The  terminology  accretion  wave  and  ero¬ 
sion  wave  refers  to  the  crest  of  the  longshore  sand 
wave  and  to  a  situation  that  has  been  observed  on 
the  California  coast  of  the  apparent  sediment 
blocking  function  of  the  crest,  whereby  a  stretch 
of  downdrift  erosion  precedes  alongshore  move¬ 
ment  of  crest  (Grove  et  al.,  1987;  Inman,  1987). 

Longshore  sand  waves  have  had  limited  study, 
despite  being  identified  early  in  the  field  of  coastal 
processes  (Evans,  1939;  Bruun,  1954;  Sonu,  1968). 
More  recently,  the  relation  of  longshore  sand 
waves  to  beachfill  response  has  heightened  interest 
in  the  study  of  this  simple  form  of  large-scale 
coastal  behavior.  Longshore  sand  waves  move  on 
a  sub-regional  or  regional  scale  and  hold  implica¬ 
tions  both  for  understanding  nearshore  morphol¬ 
ogy  change  and  for  engineering  in  the  coastal  zone, 
such  as  for  shore-protection  works  and  channel 
dredging.  Bruun  (1954)  observed  longshore  sand 
waves  along  the  Danish  North  Sea  coast  and 
distinguished  them  from  undulatory  wave  features 
in  the  shoreline  now  referred  to  as  beach  cusps. 
He  hypothesized  that  longshore  sand  waves  moved 
“in  the  direction  of  the  littoral  drift”  and  that 
inexplicable  accumulations  of  sediment  at  groin 
fields  might  be  due  to  the  arrival  of  a  longshore 
sand  wave.  Verhagen  (1989)  examined  a  100-year 
record  of  longshore  sand  waves  along  20  km  of 
Dutch  coast  and  similarly  concluded  that  periodic 
accretion  in  a  groin  field  coincided  with  the  passage 
of  longshore  sand  waves  and  not  to  classical 
trapping  of  littoral  (particulate)  drift  by  the  groins. 

Longshore  sand  waves  have  been  associated 
with  intermittency  in  sand  supply,  as  might  occur 
through  beach  nourishment  (Grove  et  al.,  1987; 
Inman,  1987),  breakup  and  movement  of  ebb-tidal 
shoals  (Bakker,  1968),  and  episodic  discharges  of 
sediment  from  rivers  (Hicks  and  Inman,  1987).  A 
longshore  sand  wave  formed  by  natural  or  artificial 
placement  of  sediment  on  the  shoreline  appears 
more  like  a  solitary  surface  water  wave  than  a 
periodic  wave,  often  having  a  crest  but  no  clear 
trough  or  wavelength.  A  series  of  longshore  sand 
waves  might  at  first  sight  look  like  a  periodic  wave, 
but  this  apparent  periodicity  is  probably  due  to  a 
cyclicity  in  sediment  supply  and  not  to  features 


arising  in  unison  and  moving  as  a  periodic  wave 
train  along  the  shoreline. 

Longshore  sand  waves  are  a  manifestation  of  the 
“collective  movement”  of  sediment  (Sonu,  1968)  in 
which  a  sand  body  maintains  morphologic  identity 
during  a  life  which  may  reach  years.  Collective 
movement  should  be  distinguished  from  movement 
of  individual  sand  particles.  Kraus  and  Horikawa 
(1990)  estimated  that  a  typical  longshore  speed  of 
a  sand  particle  is  on  the  order  of  100  times  greater 
than  the  translatory  movement  of  longshore  sand 
waves,  which  is  in  the  approximate  range  of  0.3  to 
4  km  year^‘  (Bakker,  1968;  Inman,  1987).  The 
speed  of  longshore  sand  waves  has  been  estimated 
to  be  inversely  related  to  their  length  /  as 
(Sonu,  1968,  his  fig.  11)  for  waves  with  lengths 
ranging  from  about  50  to  10,000  m.  Their  speed 
must  also  be  related  to  the  strength  and  duration 
of  the  hydrodynamic  forcing,  as  inferred  by  Evans 
(1939),  Bruun  (1954),  and  Sonu  (1968),  and  elabo¬ 
rated  on  below,  but  little  information  has  been 
available  to  study  their  movement  in  detail. 

A  quantitative  model  of  coastal  morphology 
change  must  therefore  be  capable  of  describing 
sediment  motion  at  two  scales,  relatively  rapid 
movement  of  sand  particles  and  relatively  slow 
collective  movement  of  large  sand  bodies.  As  a 
first  step  in  developing  such  an  approach,  it  is 
convenient  to  consider  longshore  sand  waves, 
which  to  a  good  approximation  undergo  2-D 
motion,  and  is  one  of  the  simplest  manifestations 
of  large-scale  coastal  behavior.  Larson  and  Kraus 
( 1991 )  showed  that  analytic  solutions  of  the  shore¬ 
line  change  equation  could  be  modified  to  account 
for  collective  movement  of  longshore  sand  waves 
by  addition  of  an  advective  term  (Inman,  1987) 
to  describe  movement  of  the  sand  form.  Earlier, 
Bakker  (1968)  examined  a  periodic  solution  of  the 
shoreline  equation  without  an  advective  term 
which  was  not  a  translatory  wave  but  described  a 
periodic  form. 

This  paper  describes  observations  and  numerical 
simulations  of  longshore  sand  waves  formed  from 
sediment  discharges  issued  from  Mecox  Inlet,  an 
ephemeral  inlet  located  on  the  southeast  shore  of 
Long  Island,  New  York,  and  connecting  Mecox 
Pond  to  the  Atlantic  Ocean.  Five  sets  of  aerial 
photographs  show  as  many  as  11  longshore  sand 
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eastern  end  of  Long  Island  which  create  longshore 
sand  waves  that  migrate  along  the  shore. 
Longshore  sand  waves  observed  updrift  of  Mecox 
Inlet  are  created  by  Sagaponack  Pond,  located 
4  km  east  of  Mecox  Bay,  intermittently  opening 
to  the  Atlantic  Ocean.  Georgica  Pond,  located 
7  km  east  of  Sagaponack  Pond,  also  releases 
migrating  sand  waves. 

The  beaches  along  the  south  shore  of  Long 
Island  are  composed  mainly  of  quartz  sand  with 
some  garnetiferous  and  magnetic  sands  and  shell 
fragments  (Leatherman  and  Allen,  1985).  The 
sand  tends  to  be  fine-to-medium  grained  with  a 
median  grain  size  of  0.25  mm  (Nersesian  et  al., 
1992).  The  beach  width  at  Southampton  ranges 
from  50  to  200  m  with  a  dune  elevation  of  approxi¬ 
mately  5  to  8  m  National  Geodetic  Vertical  Datum 
(NGVD).  Typical  beach  profiles  (Lockwood, 
Kessler  and  Bartlett,  1979)  shown  in  Fig.  2  (loca¬ 
tions  shown  in  Fig.  1),  illustrate  a  beach  face  slope 
of  approximately  1  m  vertical  to  10  m  horizontal. 
Two  bars  are  observed;  a  nearshore  bar  located 
approximately  100  m  offshore,  and  a  bar  300  to 
400  m  offshore.  Shoreline  position  was  adjusted  to 
a  common  datum  (NGVD)  by  applying  a  correc¬ 
tion  for  tide  level  at  the  time  of  photography  as 
estimated  by  the  National  Ocean  Survey  (1992) 
tide  data  and  by  translating  the  shorelines  based 
on  an  average  foreshore  slope  taken  from  the 
beach  profile  data  (Fig.  2).  The  tide  is  semi-diurnal 


with  mean  range  of  0.9  m  and  spring  range  of 
1.1  m. 

Average  waves  at  the  study  site  have  significant 
height  of  approximately  1.0  m  and  a  period  of 
about  8  sec,  according  to  a  recent  20-year  Wave 
Information  Study  (WIS)  hindcast  (Hubertz  et  al., 
1993).  The  predominant  wave  direction  is  toward 
the  west.  Seasonal  fiuctuations  include  larger 
westerly  waves  typically  occurring  from  autumn 
through  spring,  and  smaller  waves,  directed  toward 
the  east,  typically  occurring  in  summer.  The  result¬ 
ing  net  potential  longshore  transport  rate  along 
Long  Island  has  been  estimated  from  impounded 
volumes  at  jetties  to  be  300,000  m^  year”^  to  the 
west  (Panuzio,  1968).  Time  series  wave  data  used 
in  the  present  study  were  taken  from  Sep.  1956  to 
Jan.  1958  of  the  WIS  hindcast,  with  values  of  wave 
height,  period,  and  direction  available  at  3-h 
intervals.  These  data  were  selected  because  they 
were  representative  of  the  seasonal  variability  of 
wave  conditions  shown  in  the  20  years  of  WIS  data. 

3.  Observed  longshore  sand  waves 

3. 1.  Identification  of  longshore  sand  waves 

Five  recent  sets  of  vertical  aerial  photographs 
(four  sets  in  color,  one  in  black  and  white)  were 
obtained  showing  longshore  sand  waves  along 


Fig.  2.  Representative  beach  profile  shape  at  Southampton. 
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Fig.  3.  Study  reach.  Mar.  29,  1992,  showing  1 1  longshore  sand  waves. 
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Southampton  Beach  over  the  16-month  interval 
from  Sep.  1991  to  Jan.  1993.  Photographs  dated 
Sep.  4,  1991,  Dec.  20,  1991,  Dec.  18,  1992,  and 
Jan.  2,  1993  are  at  a  scale  of  1:12,000  and  were 
taken  by  TopoMetrics,  Inc.  An  additional  set  of 
photographs  taken  by  Lockwood,  Kessler,  and 
Bartlett,  Inc.,  dated  Mar.  29,  1992  is  at  a  scale  of 
1:19,200. 

Fig.  3  shows  the  locations  of  the  1 1  longshore 
sand  waves  identified  for  analysis,  giving  approxi¬ 
mately  one  wave  per  1.5  km  of  shoreline.  To  be 
considered  for  quantitative  analysis  in  this  study, 
a  longshore  sand  wave  (abbreviated  as  LSW  when 
numbered)  had  to  be  clearly  discernable  in  at  least 
three  of  the  five  sets  of  photographs.  LSW  3  was 
directly  identified  by  Terchunian  (1992)  as  created 
by  a  particular  opening  and  closing  of  Mecox 
Inlet.  It  is  noted  from  Fig.  3  that  longshore  sand 
waves  tend  to  be  asymmetric,  with  a  steeper  flank 
downdrift  of  the  predominant  direction  of  current 
(which  is  to  the  west),  similar  in  shape  to  a 


subaqueous  sand  wave  or  dune  subjected  to  uni¬ 
directional  flow.  Other  researchers  have  made  use 
of  the  analogy  between  longshore  sand  waves  and 
dunes,  including  Evans  (1939),  Bruun  (1954),  and 
Sonu  (1968).  In  the  present  case,  with  the  long¬ 
shore  current  acting  as  an  analogue  for  unidirec¬ 
tional  flow,  the  steeper  flank  lies  to  the  west. 

Fig.  4  shows  the  shoreline  and  interpreted  sub¬ 
aqueous  morphology  of  longshore  sand  waves 
inferred  from  the  aerial  photographs.  Predominant 
features  include  the  dune  line,  land-water  inter¬ 
face,  and  nearshore  morphology.  The  dune  line  is 
stationary  except  when  impacted  by  severe  storms. 
The  summer  beach  is  about  50  m  wider  than  the 
winter  beach.  The  subaqueous  morphology  associ¬ 
ated  with  sand  waves  appears  as  oblique  finger 
shoals  which  may  protrude  as  far  as  500  m  offshore 
and  are  oriented  downdrift  from  the  beach.  Similar 
finger  shoals  described  by  Sonu  (1968)  were  also 
skewed  in  the  direction  of  the  longshore  current. 
The  shoals  in  Fig.  4  are  distinguished  as  discolor- 


Fig.  4.  Plan  view  interpretation  of  nearshore  morphology. 
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ations  in  the  water  adjacent  to  many  of  the  long¬ 
shore  sand  waves  at  Southampton.  The  associated 
bathymetry  creates  a  protected  environment  for 
the  longshore  sand  waves  as  seen  by  converging 
breaking  waves  at  their  crests  (Fig.  5),  and  con¬ 
vergence  of  the  breaking  waves  is  considered  as 
being  a  major  contributor  to  the  preservation  of 
the  longshore  sand  waves. 

Three  sets  of  aerial  photographs  were  selected 
for  morphologic  analysis  and  quantification  of  the 
motion  of  the  1 1  longshore  sand  waves.  The  Sep. 
4,  1991  set  of  photographs  served  as  an  initial 
reference  condition.  Photographs  taken  on  Dec. 
20,  1991  were  chosen  to  isolate  the  winter  period 
(from  Sep.  to  Dec.)  for  determining  seasonal 
differences  in  longshore  sand  wave  migration 
speed.  The  Jan.  2,  1993  set  of  photographs  was 
selected  because  it  was  the  final  available  data  set, 
and  it  showed  less  shoreline  damage  caused  by  a 
major  storm  in  Dec.  1992  than  the  photographs 
taken  that  same  month. 

3.2.  Characterization  of  LSW geometry 

The  observed  longshore  sand  waves  were  charac¬ 
terized  as  solitary  waves.  As  shown  in  Fig.  5,  a 
wave  amplitude  rj  (maximum  distance  from  base¬ 
line  of  mean  local  shoreline  position)  and  a  wave¬ 


length  1  (distance  over  which  the  wave  protrudes 
from  the  mean  shoreline  position)  were  scaled 
from  each  of  the  three  selected  sets  of  photographs. 
An  absolute  measurement  accuracy  of  +3.5  m  can 
be  expected  from  the  technique  used.  A  threshold 
value  of  amplitude  of  20  m  was  applied  in  identi¬ 
fying  longshore  sand  waves,  representing  an  esti¬ 
mated  accuracy  in  feature  identification.  A  certain 
degree  of  interpretation  is  required  in  determining 
the  amplitude  of  the  longshore  sand  waves,  as 
some  were  occasionally  bimodal,  having  two  points 
which  could  be  identified  as  crest.  Ambiguity  in 
determining  amplitude  adds  inaccuracy  in  some  of 
the  measurements. 

A  mean  velocity  V  for  each  longshore  sand  wave 
was  determined  for  the  Dec.  20,  1991  and  Jan.  2, 
1993  data  as  the  direction  and  distance  that  the 
crest  of  the  wave  had  moved  over  the  time  elapsed 
since  Sep.  4,  1991.  For  example,  LSW  3  moved 
from  a  position  clearly  east  of  the  nearby  housing 
subdivision  on  Sep.  4,  1991  (see  reference  line  on 
Fig.  5)  to  a  location  adjacent  to  the  subdivision 
on  Dec.  20,  1991  (Fig.  6),  a  displacement  of  more 
than  200  m  in  approximately  3.5  months.  Fig.  7 
shows  LSW  3  on  Jan.  2,  1993  to  lie  west  of  the 
subdivision,  a  total  migration  distance  of  nearly 
300  m.  LSW  3  is  typical  in  that  the  mean  migration 
speed  was  three  times  greater  in  the  winter 


Fig.  5.  Longshore  Sand  Wave  3,  Sep.  4,  1991. 
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Fig.  6.  Longshore  Sand  Wave  3,  Dec.  20,  1991. 


Fig.  7.  Longshore  Sand  Wave  3,  Jan.  2,  1993. 


(Sep.  4,  1991  to  Dec.  20,  1991)  than  during  the 
16-month  period  (Sep.  4,  1991  to  Jan.  2,  1993) 
that  included  summer.  Table  1  lists  values  of  rj,  X, 
and  V  for  the  three  sets  of  photographs,  including 
mean  values  for  each  time  interval. 

LSW  10  does  not  appear  in  the  photographs  of 
Jan.  2,  1993,  and  it  is  thought  that  LSW  10  was 
overtaken  and  engulfed  by  the  faster  moving  LSW 
9.  The  photographs  indicated  that  if  a  wave  moved 


at  a  speed  0.70  km  year“^  slower  than  the  next 
consecutive  easterly  (updrift)  wave  in  the  Dec.  20, 
1991,  it  would  combine  with  the  faster  moving 
updrift  wave  to  form  one  wave  by  Jan.  2,  1993, 
similar  to  the  fate  of  LSW  10.  LSWs  5  and  8  were 
shown  in  the  Dec.  20,  1991  data  to  migrate  at 
speeds  0.75  and  0.72  km  year~^  slower  than  LSWs 
4  and  7,  respectively.  Thus,  neither  LSW  5  nor 
LSW  8  was  observed  in  the  Jan.  2,  1993  photo- 
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Table  1 

Amplitude  i/,  length  I,  and  velocity  V  of  longshore  sand  waves  observed  at  Southampton  Beach,  Long  Island,  New  York.  (Positive 
velocity  means  movement  toward  west.) 


Sep.  4, 

1991 

Dec.  20,  1991 

Jan.  2,  1993 

LSW 

r]  (m) 

X  (m) 

K(km  yr“*) 

r)  (m) 

X  (m) 

K(kmyr"') 

n  (m) 

X  (m) 

1 

56 

577 

0.37 

62 

812 

0.29 

62 

880 

2 

42 

717 

1.32 

28 

802 

0.36 

90 

1133 

3 

42 

750 

0.85 

48 

425 

0.2 

28 

558 

4 

35 

840 

1.32 

35 

111 

0.24 

77 

1165 

5 

28 

513 

0.57 

35 

398 

- 

- 

- 

6 

56 

1030 

1.18 

49 

1377 

0.43 

62 

1175 

7 

33 

527 

1.76 

49 

1358 

0.42 

62 

1422 

8 

42 

994 

1.04 

28 

1308 

- 

- 

- 

9 

42 

758 

2.18 

39 

530 

0.52 

69 

1528 

10 

35 

1217 

0.49 

56 

940 

- 

- 

- 

11 

21 

285 

0.97 

28 

342 

- 

- 

- 

Ave. 

39 

746 

1.09 

41 

779 

0.35 

64 

1130 

graphs.  LSW 

11  also  does 

not  appear  in 

.  the 

amplitude  while  decreasing  in  length,  and  three 

photographs  of  Jan.  2,  1993;  the  speed  of  LSW  11 

sand  waves  increased  in  wavelength  while  decreas¬ 

was  below  the  threshold  and  would  be  expected 

ing  in  amplitude  during  this  period.  One  sand 

to  combine  with  LSW  9  which  was  moving  more  wave  decreased  in  both  amplitude  and  length,  and 

than  1  km  year“^  faster.  LSW  1 1  was  approaching  another  sand  wave  maintained  the  same  amplitude 

the  east  jetty  at  Shinnecock  Inlet,  and  the  wave  and  decreased  in  wavelength.  These  data  suggest 

speed  was  probably  reduced  through  the  orienta-  that  no  obvious  trend  exists  between  magnitude 

tion  of  the  fillet  adjacent  to  the  jetty,  which  would  of  amplitude  and  wavelength  over  the  3.5-month 

produce  a  weaker  longshore  current  because  of  period.  However,  between  Sep.  4,  1991  and  Jan. 

smaller  angles  of  wave  breaking  to  the  local  2,  1993,  six  of  the  seven  visible  sand  waves 

shoreline.  increased  in  both  amplitude  and  wavelength.  Of 

A  relatively  strong  linear  relationship  (correla-  these  six  sand  waves,  three  would  be  expected  to 

tion  coefficient  of  0.84;  significant  to  95%  confi-  increase  in  volume  because  they  had  engulfed  one 

dence  interval)  between  observed  X  and  V  was  or  more  other  waves.  The  remaining  sand  wave, 

found  for  the  Sep.  1991-Jan.  1993  data.  These  LSW  3,  was  observed  on  Jan.  2,  1993,  to  decrease 

data  indicate  that  the  longer  longshore  sand  waves  in  both  amplitude  and  wavelength.  Overall,  the 

move  faster  than  the  shorter  ones,  in  contrast  to  16-month  observation  period  indicates  a  tendency 

the  observations  of  Sonu  (1968),  who  postulated  for  longshore  sand  waves  at  Southampton  to  main- 


that  the  longshore  velocity  of  sand  bodies  would 
vary  inversely  with  the  wavelength  as  If 

three  of  the  seven  points  denoting  combined  waves 
are  removed,  the  correlation  is  0.97  for  the  remain¬ 
ing  four  points.  No  such  trend  was  found  for  the 
Sep.  1991 -Dec.  1991,  probably  because  of  the 
short  time  period  of  observation. 

Of  the  11  longshore  sand  waves  observed  at 
Southampton,  three  increased  in  both  amplitude 
and  length  between  the  period  from  Sep.  4,  1991 
and  Dec.  20,  1991.  Three  sand  waves  increased  in 


tain  form  and  increase  in  subaerial  volume. 

Inman  (1987)  observed  that  a  longshore  sand 
wave  along  the  coast  of  California  tended  to  diffuse 
by  decreasing  in  amplitude  while  increasing  in 
length.  However,  the  longshore  sand  waves  migrat¬ 
ing  along  the  Long  Island  coast  tended  to  slightly 
increase  in  amplitude  and  length.  The  sand  wave 
observed  by  Inman  (1987)  was  produced  by  the 
artificial  injection  of  600,000  m^  of  sandy  material 
to  the  littoral  system.  Strong  dispersion  would  be 
expected  of  this  non-equilibrated  system.  In  con- 
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trast,  the  longshore  sand  waves  produced  at  Mecox 
Inlet  are,  by  the  naturally  occurring  processes, 
scaled  to  the  width  of  the  surf  zone  and  to  the 
energy  of  the  incident  waves,  suggesting  that  they 
may  be  close  to  an  equilibrium  form  that  can 
preserve  identity. 


4.  Numerical  simulation 


4.1.  Mathematical  approach  and  analytical 
formulations 


In  this  section  we  provide  background  for  the 
numerical  approach  that  follows  by  considering 
some  general  properties  of  analytical  solutions 
to  the  equation  governing  shoreline  change. 
Mathematical  modeling  of  coastal  shoreline 
change  was  originated  by  Pelnard-Considere 
(1956),  and  such  modeling  has  become  a  standard 
engineering  technique  for  predicting  the  long¬ 
term  (order  of  years  to  decades)  evolution  of  the 
position  and  shape  of  the  shoreline.  Three  basic 
assumptions  underlie  shoreline  change  models 
(Hanson  and  Kraus,  1989):  (1)  long-term  perma¬ 
nency  of  the  beach  profile  shape  (implying  the 
existence  of  an  equilibrium  form);  (2)  existence  of 
a  depth  of  closure  of  sediment  movement;  and  (3) 
dependence  of  the  sand  transport  rate  on  wave 
direction  (which  determines  the  direction  of  the 
mean  longshore  current  moving  the  sand). 

If  the  assumptions  are  made  that  the  incident 
waves  are  constant  in  height  and  direction  along¬ 
shore  and  through  time,  and  that  they  arrive  at  a 
small  angle  to  the  trend  of  the  coast,  closed-form 
analytic  solutions  are  available  (e.g.,  Pelnard- 
Considere,  1956;  Larson  et  al.,  1987).  For  this 
simple  case,  the  equation  governing  change  in  the 
shoreline  position  y,  derived  from  conservation  of 
sand,  reduces  to  the  1-D  diffusion  equation: 


dt  ^  dx^ 


(1) 


in  which  the  shoreline  position  y  is  a  function  of 
the  distance  x  along  the  shore  and  the  time  t,  and 
e  is  an  empirical  coefficient  related  to  the  particular 
predictive  longshore  sand  transport  rate  formula 


and  the  depth  of  closure  (depth  to  which  the 
beach  profile  moves  when  the  shoreline  moves). 
Numerous  useful  closed-form  or  analytical  solu¬ 
tions  can  be  found  for  Eq.  1  (Larson  et  al.,  1987); 
however,  this  equation  describes  diffusion  which 
tends  to  obliterate  the  distinct  and  persistent  shore¬ 
line  sand  forms  which  are  addressed  in  this  study. 

In  recent  years,  interest  has  heightened  in  what 
Sonu  (1968)  called  collective  movement  of  sand, 
in  which  large  bodies  having  the  dimensions  on 
the  order  of  the  width  of  the  surf  zone  move 
alongshore.  The  longshore  movement  of  large  sand 
forms  can  be  incorporated  in  Eq.  1  by  including  a 
form  advective  term  V{dy/dx)  (Inman,  1987; 
Larson  and  Kraus,  1991),  which  can  be  derived 
as  in  forming  the  standard  advection-diffusion 
equation  for  a  conservative  substance,  to  give: 


dt^^dx-^dx^ 


(2) 


in  which  V  is  the  migration  speed  of  the  sand  wave 
and  is  presumably  a  function  of  the  hydrodynamic 
conditions,  which  vary  in  time  and  also  depend  on 
the  offshore  bathymetry. 

No  predictive  expressions  are  available  to  specify 
V,  although  Sonu  (1968)  deduced  that  the  migra¬ 
tion  speed  of  “sinuous  shorelines  and  crescentic 
bars”  was  a  function  of  spacing  of  the  forms,  with 
larger  bodies  moving  more  slowly.  Sand  wave 
speed  on  the  southern  coast  of  California  have 
been  reported  in  the  range  of  0.5  to  4  km  year“^ 
(Grove  et  al.,  1987;  Inman,  1987).  Similar  values 
have  been  found  by  others  for  large  sand  protru¬ 
sions  moving  alongshore  and  establish  the  order 
of  magnitude  for  V.  The  speed  V  is  taken  as  a 
constant  for  obtaining  simple  analytic  solutions 
for  Eq.  2.  As  shown  by  Larson  and  Kraus  (1991), 
if  the  transformation  is  made  to  a  coordinate 
system  moving  in  the  positive  x  direction  with 
speed  Fby  defining  arguments  <j>=Vt—x  and  0  = 
t,  then  by  the  chain  rule  of  differentiation,  Eq.  2 
reduces  to  the  same  form  as  Eq.  1  (let  ^  go  to  x 
and  6  go  to  t).  Therefore,  if  there  are  no  distur¬ 
bances  at  the  boundaries,  analytical  solutions  of 
Eq.  1  are  also  solutions  of  Eq.  2,  but  with  the 
arguments  ^  and  6.  Eq.  2  thus  describes  the  migra¬ 
tion  of  a  sand  body  with  a  characteristic  velocity 
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V  that  tends  to  diffuse  and  lose  identity.  A  numeri¬ 
cal  solution  is  required  to  achieve  more  generality 
such  as  feedback  of  bathymetry  to  the  forcing 
waves  and  to  allow  time-varying  and  longshore- 
varying  wave  conditions,  described  in  the  next 
section. 

In  Eq.  2,  the  term  proportional  to  V  represents 
the  movement  of  a  shoreline  form,  such  as  long¬ 
shore  sand  waves,  and  the  term  on  the  right  side 
represents  shoreline  change  produced  by  particu¬ 
late  motion.  Let  us  set  £  =  0,  implying  no  individual 
particulate  motion,  giving: 


8y  By 

—  +  V  — 
Bt  Bx 


=  0 


(3) 


A  periodic  or  wave-like  solution  (longshore  sand 
wave)  to  this  equation  is: 

/  2n  2n 

>'  =  '?maxCOSl  y  X— 

in  which  is  the  amplitude  of  the  sand  wave, 
X  is  the  wavelength  of  the  sand  body,  and  t  is  the 
period  of  the  wave. 

Closed-form  solutions  require  further  simplifi¬ 
cations,  the  major  ones  being  a  schematic  initial 
shape  of  the  shoreline  feature  and  waves  of  con¬ 
stant  height  and  direction  through  time  and  along¬ 
shore.  Therefore,  to  allow  generality  and  to  more 
realistically  represent  time-varying  waves,  wave 
transformation  across  shore,  and  migration  and 
evolution  of  a  longshore  sand  wave,  a  numerical 
shoreline  change  simulation  model  was  developed. 


4.2.  Shoreline  change  and  sand  transport 


numerical  model: 


By  By 

—+  V  — 
Bt  Bx 


/)<•  Bx 


(5) 


The  second  term  on  the  left  side  of  Eq.  5  is  an 
advection  term  that  governs  the  movement  of  the 
longshore  sand  wave. 

The  dependence  of  V  on  possible  governing 
variables  has  not  been  previously  examined,  and 
an  estimation  method  to  calculate  it  in  the  model 
had  to  be  developed  in  the  present  study.  One 
expects  with  Bruun  (1954)  and  others  that  a 
longshore  sand  wave  should  move  in  the  direction 
of  the  littoral  drift  or  with  the  longshore  current. 
One  would  also  expect  that  a  stronger  current 
would  cause  more  rapid  movement.  What  is 
“strong?” 

In  measurements  of  sand  waves  in  uniform  flow 
of  the  upper  part  of  an  estuary,  Nasner  (1974,  his 
fig.  5)  found  a  well-defined  linear  relation  between 
the  speed  of  sand  waves  (or  dunes)  in  the  Weser 
River,  Germany,  and  the  mean  freshwater  dis¬ 
charge  that  occurred  between  depth  soundings. 
Kraus  and  Dean  (1987)  found  correlation  between 
the  measured  particulate  longshore  sand  transport 
rate  and  the  longshore  discharge  of  water  in  the 
surf  zone  at  Duck,  North  Carolina,  facing  the 
Atlantic  Ocean.  Here,  we  associate  the  speed  of  a 
longshore  sand  wave  with  the  longshore  discharge 
of  water  rather  than  with  longshore  current  speed 
(such  as  the  mean  current  speed). 

Following  Kraus  and  Dean  (1987),  a  longshore 
water  discharge  parameter  R  is  defined  as: 


If  the  particulate  transport  rate  is  given  by 
the  “CERC”  formula  (Shore  Protection  Manual, 
1984),  then  Qp  =  Qo  sm(2ff^,)  in  which  is  the 
angle  of  the  breaking  wave  crests  to  the  trend  of 
the  local  shoreline,  and  Qg  the  amplitude  of  the 
longshore  sand  transport  rate  (Larson  et  al.,  1987). 
With  this  choice  of  the  longshore  sand  transport 
rate  formula,  e  becomes  e  =  2QglD^  in  which 
Dc  is  the  depth  of  closure,  and  d^yldx^  = 
{2Qg)~^BQpldx.  The  preceding  substitutions  into 
Eq.  2  result  in  the  equation  governing  the 
movement  of  longshore  sand  waves  used  in  the 


^  =  2^bJ'b'’is  (6) 

in  which  d^  is  the  depth  at  wave  breaking,  yf,  is 
the  distance  from  the  shoreline  to  the  break  point, 
and  v,s  is  the  mean  velocity  of  the  longshore 
current.  The  wave  transformation  model  takes 
waves  over  locally  straight  and  parallel  contours 
from  offshore  to  the  break  point,  and  the  assump¬ 
tion  of  a  beach  profile  shape  allows  calculation  of 
and  jb-  The  longshore  current  velocity  can  also 
be  calculated  from  parameters  given  by  the  wave 
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model  as; 

=  y-J^b  sin(2eb)  (7) 

Eq.  7  is  an  empirical  expression  found  by  Komar 
and  Inman  (1970)  in  their  longshore  sand  trans¬ 
port  tracer  field  experiments,  in  which  y  is  ratio  of 
wave  height  to  water  depth  at  breaking,  g  is 
acceleration  due  to  gravity,  and  6^  is  angle  of  wave 
crests  at  breaking  to  the  local  shoreline. 

The  volume  rate  of  transport  Qlsw  of  a  long¬ 
shore  sand  wave  can  be  estimated  from  the  data 
set  as: 

(SlSW  )est.m  =  ^7)  F  (8) 

in  which  rj  is  the  amplitude  of  the  sand  wave  at  a 
given  calculation  time  step,  and  the  quantity  D  is 
the  sum  of  the  depth  of  profile  closure  and  the 
elevation  of  the  active  berm.  The  closure  depth 
was  taken  as  6  m,  and  the  elevation  of  the  berm 
was  2  m. 

The  relation  between  the  volume  rate  of  long¬ 
shore  sand  wave  migration  and  the  longshore 
discharge  is  taken  as: 

(2LSw)calc  =  ®(7^~7?<;ri,)  (9) 

and  defines  an  empirical  proportionality  coefficient 
a.  An  estimate  for  the  empirical  coefficient  a  was 
obtained  by  substituting  Eq.  8  into  Eq.  9  and  using 
representative  values  for  observed  sand  wave  speed 
and  the  hindcast  wave  conditions,  giving  a  = 
8.8  X  lO^"*.  This  value  represents  an  average  over 
summer  and  winter  wave  conditions.  The  threshold 
value  7?crit  is  unknown;  for  simplicity  it  was  set 
equal  to  2.4  m^  s“*  found  for  particulate  transport 
( Kraus  and  Dean,  1987).  for  collective  motion 
is  probably  larger,  based  on  the  results  of  Nasner 
(1974). 

In  the  model,  V  is  calculated  at  each  time  step 
from  calculated  breaking  wave  parameters  and  the 
amplitude  of  the  sand  wave  at  that  time  step, 
allowing  the  wave  conditions  to  control  the  speed 
and  direction  of  movement  of  the  longshore  sand 
wave  on  a  time-dependent  basis. 


4.3.  Wave  model 

The  wave  model  of  Kraus  and  Harikai  (1983) 
was  employed.  At  each  calculation  time  step,  the 
model  refracts  and  shoals  waves  of  specified  height 
and  direction  over  locally  straight  and  parallel 
contours  to  breaking  using  the  breaking  criterion 
H^  =  0.7S  dy,,  in  which  is  the  wave  height  at 
breaking.  An  important  feature  of  the  wave  model 
is  what  Kraus  and  Harikai  called  a  “contour 
correction”  to  Snell’s  Law.  By  means  of  the  con¬ 
tour  correction,  the  orientation  of  the  shoreline  at 
each  grid  point  alongshore  is  introduced  into 
Snell’s  Law  to  the  depth  of  closure;  in  the  present 
situation,  the  contour  correction  in  effect  replicates 
the  curvature  of  the  subaerial  portion  of  the  sand 
wave  in  the  offshore  depth  contours  out  to  deep 
water  where  the  waves  are  no  longer  effected  by 
the  bathymetry.  As  shown  in  Fig.  4,  the  bathyme¬ 
try  associated  with  the  longshore  sand  waves  is 
more  complex  than  the  simple  assumption  involv¬ 
ing  the  shoreline  shape,  but  the  contour  correction 
at  least  partially  accounts  for  feedback  between 
the  local  subaqueous  morphology  and  the  waves. 

The  contour  correction  causes  waves  to  refract 
toward  a  longshore  sand  wave,  much  as  waves 
tend  to  converge  toward  a  headland.  Convergence 
reduces  the  relative  angle  between  the  breaking 
waves  and  local  shoreline,  and  as  shown  in 
Section  5,  significantly  reduces  the  degree  of  sand 
wave  dispersion  in  the  model  runs. 

4. 4.  Numerical  solution  scheme 

Eq.  5  was  solved  by  an  explicit  finite-difference 
method.  Considerable  effort  was  spent  in  determin¬ 
ing  the  best  representation  of  the  advection  term 
describing  movement  of  the  longshore  sand  wave. 
In  independent  tests,  this  term  was  represented  by 
central  differencing,  linear  upwind  differencing, 
and  the  quick  and  quickest  (quadratic)  upwind 
difference  methods  (Leonard,  1979).  Without  par¬ 
ticulate  longshore  transport  (2?  =  0),  allowing  pure 
advection  of  a  longshore  sand  wave,  the  central 
difference  scheme  showed  persistent  small  oscilla¬ 
tions  at  the  base  of  the  calculated  wave,  where  it 
joined  to  straight  shoreline,  and  standard  linear 
upwind  differencing  introduced  severe  numerical 
diffusion.  Of  the  quick  and  quickest  methods,  the 
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Distance  Alongshore,  km 


Fig.  8.  Calculated  longshore  sand  wave  evolution  without  contour  correction.  (Full  equation,  no  contour  correction,  H=  \  m, 
r=  8  s,  direction  20°.) 


QUICKEST  method  proved  most  stable  and  had  little 
numerical  diffusion.  For  pure  advection  of  a  long¬ 
shore  sand  wave  the  size  of  LSW  3  under  constant 
obliquely  incident  waves  (1  m,  8  s,  20°  angle  at 
20-m  depth)  for  4  months,  the  crest  of  the  long¬ 
shore  sand  wave  decreased  about  0.5  m.  In  the 
tests  a  time  step  of  3  h  and  uniform  grid  spacing 
of  50  m  were  used.  Based  on  the  test  results, 
QUICKEST  was  employed. 

5.  Results 

In  reproducing  observed  longshore  sand  wave 
movement,  a  time  step  of  1  h  and  a  grid  spacing 
of  50  m  were  used.  At  each  time  step,  a  check  was 
made  to  see  if  the  Courant  number  exceeded  0.5 
as  described  by  Kraus  and  Harikai  (1983).  If  the 


value  0.5  was  exceeded,  the  time  step  was  reduced 
by  factors  of  1/2  until  the  condition  was  satisfied. 
In  this  manner,  arbitrary  incident  wave  conditions 
as  given  by  the  WIS  hindcast  could  be  used  with 
arbitrary  shoreline  configurations  in  automated 
manner  without  producing  numerical  instability. 

In  this  study,  the  movement  of  a  single  longshore 
sand  wave  was  investigated  for  clarity.  The  length 
of  shoreline  was  taken  as  18  km  to  eliminate 
boundary  effects,  and  sand  was  allowed  to  move 
on  and  off  of  the  grid  by  particulate  transport. 
The  dimensions  of  the  initial  form  of  the  sand 
wave  were  that  of  LSW  3.  The  position  of  the 
shoreline  was  set  at  100  m  from  an  arbitrary 
baseline  for  ease  of  reading  the  plots.  In  the  figures 
to  follow,  the  reader  should  be  aware  of  the  large 
vertical  exaggeration  (approximately  350  to  1  in 
most  figures). 
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Distance  Alongshore,  km 


Fig.  9.  Calculated  longshore  sand  wave  evolution  with  contour  correction.  (Full  equation,  H=  1  m,  7=8  s,  direction  20°.) 


5.1.  Model  tests 

Many  tests  were  performed  to  verify  model 
predictions  under  symmetry  of  incident  wave 
action,  positive  longshore  sand  wave  shape,  and 
reversed  or  negative  wave  shape  (pure  erosion 
wave).  The  model  performed  correctly  in  these 
tests. 

Fig.  8  shows  calculated  longshore  sand  wave 
migration  under  constant  incident  wave  action 
with  the  wave  contour  correction  turned  off.  The 
incident  waves  in  10-m  depth  had  significant  height 
of  1  m,  period  of  8  s,  and  direction  of  20°.  With 
these  strong  and  steady  waves,  after  two  years  (an 
idealized  year  of  360  days  was  used)  the  crest  of 
the  longshore  sand  wave  moved  about  6  km.  The 
entire  sand  wave  body  moved  west,  but  there  was 
significant  decrease  in  amplitude,  from  34.2  to  8.3 
m,  and  extraordinary  spreading  was  produced  by 


the  particulate  transport  term.  The  substantial 
spreading  and  decrease  in  amplitude  are  not  in 
accord  with  the  field  observations. 

The  same  case  as  above  was  run  with  the  wave 
contour  correction  turned  on,  and  the  results  are 
shown  in  Fig.  9.  Compared  to  the  case  with  the 
contour  correction  inoperative,  the  body  of  long¬ 
shore  sand  wave  remains  much  more  intact  and 
its  amplitude  decreases  less,  from  34.2  to  16.5  m. 
In  addition,  the  longshore  sand  wave  becomes 
asymmetrical,  possessing  a  steeper  flank  downdrift 
in  qualitative  agreement  with  the  field  observa¬ 
tions.  Based  on  these  and  similar  results,  the 
contour  correction  is  employed  in  subsequent 
calculations. 

The  next  test  shown  is  that  for  two  half-year 
seasons,  with  summer  waves  for  the  first  six  months 
of  height  0.5  m,  period  8  sec,  and  direction  —20° 
(waves  directed  to  the  east)  and  winter  waves  of 
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height  1  m,  period  8  s,  and  direction  20°  (waves 
directed  to  the  west).  The  calculated  results,  output 
at  4-month  intervals,  are  shown  in  Fig.  10.  For  the 
first  four  months  of  pure  mild  summer  waves,  the 
longshore  sand  wave  moves  less  than  0.5  km  to 
the  east,  and  its  decrease  in  amplitude  is  relatively 
small.  The  next  4-month  interval  includes  two 
months  of  summer  waves  followed  by  two  months 
of  winter  waves.  The  sand  waves  therefore  experi¬ 
ences  a  net  movement  to  the  west,  and  the  ampli¬ 
tude  decreases  somewhat  more  than  in  the  first 
four  months.  The  final  4-month  increment  contains 
purely  winter  waves,  and  the  crest  of  the  sand 
wave  moves  the  most,  approximately  1.3  km.  This 
example  shows  that  sand  waves  can  move  either 
east  or  west  in  the  model,  and  that  the  translation 
speed  and  spreading  of  the  sand  wave  are  related 
to  the  incident  wave  characteristics. 


5.2.  Simulation  of  observed  sand-wave  movement 

As  previously  mentioned,  comparisons  of  calcu¬ 
lated  and  observed  evolution  of  only  one  longshore 
sand  wave  are  presented  for  clarity.  The  period 
Sep.  1956-Jan.  1958  was  used  as  input  to  the 
model,  and  incident  waves  were  input  for  the 
16-month  observation  period.  As  shown  in  Fig.  11, 
the  observed  sand  wave  (LSW  3)  did  not  decrease 
in  amplitude  for  the  first  four  months,  and  it 
appeared  to  become  narrower.  For  the  same 
increment,  the  calculated  sand  wave  moved  in  the 
correct  direction  (as  it  did  for  all  time  increments), 
but  not  quite  as  far.  In  addition,  the  calculated 
wave  spread.  The  observed  sand  wave  on  Jan.  2, 
1993  had  a  double  peak  at  its  crest,  and  the  crest 
of  the  final  calculated  wave  lies  in  the  approximate 
center  of  these  two  peaks.  Overall,  reasonable 


Fig.  10.  Calculated  longshore  sand  wave  evolution  for  idealized  seasonal  waves.  (Artificial  seasons,  Simi=0.5  m,  8  s,  —20°; 
Win  =  l  m,  8  s,  20°.) 
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quantitative  agreement  is  seen,  but  with  the  calcu¬ 
lated  sand  wave  being  too  broad. 

As  shown  in  Fig.  11,  the  observed  sand  waves 
were  much  sharper  in  shape  (seeming  unnatural 
due  to  the  vertical  exaggeration  of  the  plots)  than 
the  calculated  waves.  The  initial  form  for  the 
calculated  sand  wave  was  prepared  from  a  nar¬ 
rower  wave  that  was  subjected  to  incident  wave 
action  to  obtain  the  correct  initial  amplitude;  the 
extra  width  and  smoothness  result  from  the  spread¬ 
ing  by  particulate  transport.  A  general  characteris¬ 
tic  of  the  model  was  to  produce  broader  forms 
than  observed. 


6.  Concluding  discussion 

In  this  study,  multiple  longshore  sand  waves 
were  clearly  tracked  along  Southampton  Beach. 


The  sand  waves  had  an  average  length  of  0.75  km 
and  average  amplitude  of  40  m.  Migration  speed 
was  greatest  under  winter  waves  and  reached  an 
average  1.09  km  year~^  in  winter,  whereas  yearly 
average  speed  including  summer  calms  was  a  much 
lower  0.35  km  year“^  During  a  16-month  obser¬ 
vation  period,  the  longshore  sand  waves  did  not 
disperse  during  steady  migration  westward  in  the 
direction  of  predominant  wave  incidence,  and 
longer  waves  tended  to  move  faster  than  those 
with  smaller  wavelengths. 

A  numerical  simulation  model  was  developed 
that  reproduced  many  of  the  dynamic  properties 
of  longshore  sand  waves.  The  model  contains 
transport  contributions  describing  both  individual 
(particulate)  grain  movement  and  collective  move¬ 
ment  of  the  sand  body.  The  speed  and  direction 
of  the  sand  waves  were  related  to  the  longshore 
discharge  of  water,  thereby  giving  the  speed  in  a 
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predictive  equation  in  terms  of  the  incident  wave 
conditions,  sand  wave  amplitude  and  depth  of 
closure.  The  single  empirical  coefficient  in  the 
transport  rate  formula  was  determined  based  on 
measured  and  calculated  average  conditions.  The 
calculated  sand  waves  moved  more  slowly 
in  summer  than  winter,  in  agreement  with 
observations. 

An  essential  feature  of  the  model  is  use  of  a 
contour  correction,  whereby  the  depth  contours  of 
the  longshore  sand  wave  are  extended  offshore, 
entering  the  wave  refraction  calculation.  The  con¬ 
tour  correction  causes  incident  water  waves  to 
converge  on  the  sand  wave,  thereby  reducing  the 
sand-wave  dispersion.  The  contour  correction  also 
reproduced  the  observed  asymmetry  in  longshore 
sand  wave  form,  which  tends  to  show  a  steeper 
downdrift  flank. 

Refraction  may  also  explain  the  longevity  of 
observed  sand  waves.  However,  the  calculated 
sand  wave  evolution  showed  more  decrease  in 
amplitude  and  more  spreading  than  the  observed 
sand  waves.  It  is  believed  that  one  possible  reason 
for  less  spreading  in  nature  is  that  the  subaqueous 
portion  of  the  sand  wave  is  more  extensive  than 
inferred  from  the  subaerial  portion.  If  this  is  true, 
than  the  basic  shoreline  modeling  approach  will 
have  to  be  revised  to  allow  non-uniformity  in 
beach  contours.  Bathymetric  surveys  of  longshore 
sand  waves  need  to  be  made  to  provide  informa¬ 
tion  for  improving  our  quantitative  understanding 
and  predictive  capability. 

The  present  study  indicates  that  successful 
numerical  simulation  of  longshore  sand  waves  is 
likely  possible,  but  detailed  field  surveys  are 
required  to  provide  needed  information  to  progress 
in  the  mathematical  description.  Longshore  sand 
waves  are  worthy  of  additional  study  because  they 
are  a  simple  form  of  collective  sand  body  motion 
and  because  they  are  of  engineering  significance  to 
such  projects  as  shore  protection  and  dredging  of 
navigation  channels. 
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